Thermal radiation from an expanding viscous medium 
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The effects of viscosity on the space time evolution of QGP produced in nuclear collisions at RHIC 
energies have been studied. The entropy generated due to the viscous motion of the fluid has been 
taken into account in constraining the initial temperature by the final multiplicity (measured at the 
freeze-out point). The viscous effects on the photon spectra has been introduced consistently through 
the evolution dynamics and phase space factors of the particles participating in the production 
process. We notice a stronger effect on the photon spectra originating from QGP than hadronic 
matter. A detectable shift is observed in the space-time integrated pr distribution of photons due 
to dissipative effects. 
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I. INTRODUCTION 

Nuclear collisions at Relativistic Heavy Ion Col- 
lider (RHIC) and Large Hadron Collider (LHC) en- 
ergies are aimed at creating a thermalized state 
of quarks and gluons called quark gluon plasma 
(QGP). The weakly interacting picture of the QGP 
stems from the perception of asymptotic freedom of 
QCD at high temperatures and densities. However, 
the experimental data from RHIC [l| , especially the 
measured elliptic flow of hadrons indicate that the 
matter produced in Au+Au collisions exhibit prop- 
erties which are more like a strongly interacting 
liquid than a weakly interacting gas. The magni- 
tude of the transport coefficients can be used to un- 
derstand the strength of the interaction within the 
QGP. Therefore, the study of the transport proper- 
ties of QGP and hot hadrons is of paramount impor- 
tance in characterizing the matter formed in heavy 
ion collisions (HIC) at relativistic energies. For ex- 
ample, the shear viscosity or the internal friction 
of the fluid symbolizes the ability to transfer mo- 
mentum over a distance of about a mean free path. 
Therefore, in a system where the constituents inter- 
act strongly the transfer of momentum is performed 
easily - resulting in lower values of r\. Consequently 
such a system may be characterized by a small value 
of rj/s where s is the entropy density. On the other 
hand, for a weakly interacting system the momen- 
tum transfer between the constituents become stren- 
uous which gives rise to large 77. The importance of 
viscosity also lies in the fact that it damps out the 
variation in the velocity and makes the fluid flow 
laminar. A very small viscosity (large Reynold num- 
ber) may make the flow turbulent. A lower bound on 
the value of -q/s has been found using AdS/CFT Q 
(see also Q). 

Collisions between nuclei at ultra-relativistic en- 
ergies produce charged particles - either in the 
hadronic or in the partonic state, depending on the 
collision energy. Interactions among these charged 



particles produce real photons. Because of their 
nature of interaction, the mean free path of pho- 
tons in the medium (hadronic or partonic) is large 
compared to the size of the system formed in HIC. 
Therefore, photons emanating from such a system 
brings out the information of the source point very 
efficiently 4-6>] (see 0-Q for review) and hence elec- 
tromagnetic probes (photons and lepton pairs) may 
play crucial role in extracting the transport coeffi- 
cients. 

The effects of viscosity on the photon spectra re- 
sulting from HIC enter through two main factors: 
(i) the modification of the phase space factor due 
to the deviation of the system from equilibrium and 
(ii) the space time evolution of the matter governed 
by dissipative hydrodynamics. One more important 
issue deserves to be mentioned here. Normally, the 
initial temperature (Tj) and the thermalization time 
(Ti) are constrained by the measured hadron multi- 
plicity (dN/dy) . This approach is valid for a system 
where there is no viscous loss and the time reversal 
symmetry is valid. However, for a viscous system 
the entropy at the freeze-out point (which is propor- 
tional to the multiplicity) contains the initially pro- 
duced entropy as well as the entropy produced dur- 
ing the space time evolution due to non-zero shear 
and bulk viscosity. Therefore, the amount of entropy 
generated during the evolution has to be subtracted 
from the total entropy at the freeze-out point and 
the remaining part which is produced initially should 
be used to estimate the initial temperature. There- 
fore, for a given dN/dy (which is associated with the 
freeze-out point) and t% the magnitude of Ti will be 
lower in case of viscous dynamics compared to ideal 
flow. 

Effects of viscosity on the transverse momen- 
tum distribution of photons was earlier considered 
in [10|, [ill and recently the interest in this field is 
renewed p^ - fl4| . Beyond a certain threshold in col- 
lision energy the system is expected to be formed in 
QGP phase which will inevitably make a transition 



to the hadronic matter later. The measured spectra 
contain contributions from both QGP and hadronic 
phases. Therefore, it becomes imperative to esti- 
mate the photon emission with viscous effects from 
QGP as well as hadrons and identify a kinematic 
window where photons from QGP dominate. While 
in some of the earlier works [L2H14J contributions 
from hadrons were ignored, in others [10|, UM the ef- 
fects of dissipation on the phase space factors were 
omitted. In the present work we study the effects 
of viscosity on the thermal photon spectra originat- 
ing from QGP and hadronic matter and argue that 
photons can be used as a very useful tool to estimate 
r\l s and hence characterize the matter. 

The paper is planned as follows. In the next sec- 
tion the kinetic theory formalism for evaluating the 
photon emission rate in the independent particle ap- 
proximation [l5| is discussed. In section III we dis- 
cuss the viscous effects on the phase space distri- 
butions of the partons or hadrons participating in 
the photon production processes. Hydrodynamical 
evolution with viscous effects has been mentioned in 
section IV. Section V is devoted to results and sec- 
tion VI is dedicated to the summary and discussion. 



II. PRODUCTION OF THERMAL 
PHOTONS 

The transverse momentum (pt) distribution of 
photons from a reaction of the type: 1 + 2 — > 3 + 7 
taking place in a thermal bath at a temperature, T 
is given by [15j : 
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where R is the rate of photon production per unit 
four- volume, M is the over all degeneracy for the re- 
action under consideration, p,, Ei and fi(Ei) are the 
three momentum, energy and thermal phase space 
factor of the particle i (either parton or hadron). 
\M\ 2 is the square of the invariant amplitude for 
the process under consideration. After some straight 
forward algebra Eq. Q] can be simplified to (see Ap- 
pendix A): 
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The collision of nuclei at RHIC and LHC energies 
is expected to produce QGP. This state of matter, 
once created with high internal pressure will undergo 
rapid expansion and consequently will cool down to 



the temperature, T c for QGP to hadron transition. 
In a first order phase transition scenario the sys- 
tem remains in a mixed phase of QGP and hadrons 
at T c until a time where the entire QGP converts 
to hadrons. The thermal equilibrium may also be 
maintained in the hot hadronic phase till the freeze- 
out point (achieved at a temperature, Tp) where the 
mean free path of the hadrons is too large for colli- 
sions to take place. 

The measured photon spectra (dN/d 2 pTdy) is the 
yield obtained after performing the space time inte- 
gration over the entire evolution history - from the 
initial state to the freeze-out point. Therefore, Eq.[5] 
needs to be integrated over the four volume to con- 
nect the theoretical results with experiments: 
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where i = Q,M,H represents QGP, mixed (coex- 
isting phase of QGP and hadrons) and hadronic 
phases respectively. The effects of viscosity enter 
the photon spectra through the space time evolu- 
tion governed by the dissipative hydrodynamics and 
the phase space factor, /j's in Eq. [2] 



A. Thermal photons from QGP 

The contribution from QGP to the spectrum of 
thermal photons due to annihilation {qq—^gj) and 
Compton (q(q)g -» q(q)j) processes has been calcu- 
lated in [la 11711 using hard thermal loop (HTL) ap- 
proximation [18J. Later, it was shown that photons 
from the processes [l9j: gq^-gq^f, qq^-qqj, qqq^-qi 
and gqq^f-gj contribute in the same order 0(aa s ) 
as the Compton and annihilation processes. The 
complete calculation of emission rate from QGP to 
order a s has been performed by resumming ladder 
diagrams in the effective theory [201 ]. However, in 
the present work we consider only the Compton and 
annihilation processes for photon production. We 
expect that the shift in the photon spectra from the 
ideal to the viscous scenario will not alter drastically 
with the replacement of the Compton 4- annihilation 
rates by the rate obtained in Ref. |20| . 



B. Thermal photons from hadronic matter 

A set of hadronic reactions with all possible iso- 
spin combinations have been considered for the pro- 
duction of photons |2l| - [24| from hadronic matter. 
The relevant reactions and decays for photon pro- 
duction are: (i) tttt — > pj, (ii) n p — ► nj (with 
7r, p, w, (f> and a\ in the intermediate state J23|). 
(iii)7T7r — > rjj and (iv) iri] — > ttj, p — >■ irirj 



and w — >■ it"/. The corresponding vertices's are ob- 
tained from various phenomenological Lagrangians 
described in detail in Ref. [2lH23j. The effect of 
hadronic dipolc form factors has been taken into ac- 
count in the present work as in |24| . 



III. VISCOUS CORRECTION TO THE 
DISTRIBUTION FUNCTION 
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where p' z — rriT sinh(y — rf) is the z-componcnt of the 
momentum in the fluid co-moving frame. The phase 
space distribution with viscous correction (|5J| thus 
enters the production rate of photon through Eq. [2] 



We assume that the system is slightly away 
from equilibrium which relaxes back to equilibrium 
through dissipative processes. Here we briefly recall 
the main considerations leading to the commonly 
used form for the first viscous correction, <5f to the 
phase space factor, / defined as follows [25j : 
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where /o is the equilibrium distribution function, 
(V a up) = V a U[} + Vf3U a - |A Q/3 V 7 -u 7 , A Q 4 
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g a p - u a up, V Q = (g a p - UaU^d 13 , m m being the 
four-velocity of the fluid. The coefficients C and A 
can be determined in the following way. Substitut- 
ing / in the expression for stress-energy tensor T^ v 
we get, 
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where T{f v = (e + P)u^u u — g^ v P is the energy mo- 
mentum tensor for ideal fluid. From general con- 
siderations |26| the dissipative part can be written 
as 
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Equating the part containing Sf from (g|) with ^ , C 
and A can be expressed in terms of the coefficients 
of shear (77) and bulk (£) viscosity respectively in 
terms of which the phase space distribution for the 
system can be written as: 
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For a boost invariant expansion in (1+1) dimen- 
sion this can be simplified to get, 



where 
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IV. EXPANSION DYNAMICS 

As mentioned before the px distribution of ther- 
mal photons is obtained by integrating the emis- 
sion rate over the evolution history of the expand- 
ing fluid. Relativistic viscous hydrodynamics can be 
used as a tool for the space-time dynamics of the 
fluid. 




x(fm) 

FIG. 1: Variation of temperature with proper time for 
different phases for various values of the shear viscosities. 
Inset shows the effect of viscosity on the cooling of the 
QGP phase (in an amplified scale) for different values of 
T)/s. 



For a (1+1) dimensional boost invariant expan- 
sion [27| the evolution equation d^T^ = 0, can be 

written as 1281: 



de 



— = b + C)/T 
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where P is the pressure and e is the energy density. 
We assume that the baryonic chemical potential is 
small in the central rapidity region for RHIC/LHC 
collision energies. Therefore, the equation corre- 
sponding to the net baryon number conservation 
need not be considered in these situations. 

We assume that the system achieves thermal equi- 
librium at a time r» after the collision at an ini- 
tial temperature T. With this initial condition and 



equation of state (EoS) P 
can be written as flOf : 



e/3 the solution of Eq. |H] 
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where ^q=(|?/qo + Cqo), r)Q0=riQ/ T3 = 4 aQ(v/s)Q 
andC Q o-CQ/r 3 =4a Q (C/ S )Q. 

Eq.[TU]dictates the cooling of the QGP phase from 
its initial state to the transition temperature, T c at 
a time, r 9 , when the QGP phase ends. 

In a first order phase transition scenario, the pure 
QGP phase is followed by a coexistence phase of 
QGP and hadrons. The energy density, shear and 
bulk viscosities in the mixed phase can be written in 
terms of the corresponding quantities of the quark 
and hadronic phases at temperature T c as follows 
.10]: 

cm(t) = f Q e Q (T c ) + (l-f Q (r))e H (T c ) 
Vm(t) = fQVQ(Tc) + (l-f Q (T))m(Tc) 
Cm(t) = IqCq{T c ) + (1-I q {t))Ch{T c ) (11) 

where Jq{t) (/h(t)) indicates the fraction of the 
quark (hadronic) matter in the mixed phase at 
a proper time r. We have €Q(T c )=3aQT^ + B, 
e H (T c )=3a H T?, a Q = g Q ir 2 /9Q, a H = 5ff7r 2 /90, 
B is the bag constant, gQ (gn) denote statistical 
degeneracy for the QGP (hadronic) phase. In the 
mixed phase the temperature remains constant but 
the energy density varies with time as the conversion 
of QGP to hadrons continues. This time variation 
is executed through Jq(t). Substituting Eqs. [TT1 in 
Eq. [5] and solving for /q(t) we get [1C| . 



Iq = 



,-b/T 



TQ 



CC 



b/r' 



ob/r' 



dT ' + lQ e (b/r Q -b/r) 



where a=4e ff /(3Ae), b=[4(?7Q 

i 4 



(12) 
tih)/3 + 2(C Q - 
Cff)]/Ae, c= (|tjh + 2C ff )/Ae and Ae=e Q - e H . 
Eq. [T2] indicates how the fraction of QGP in the 
co-existence phase evolves with time. 

The variation of T with r in the hadronic phase 
can be obtained by solving Eq. |H] with the bound- 
ary condition T — T c and r — th, where th is the 
(proper) time at which the mixed phase ends i.e. 
when the conversion of QGP to hadronic matter is 
completed, 
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Similar to QGP, P = e/3 has been used for hadronic 
phase. For a vanishing bulk viscosity (£ = 0) the 
cooling of the QGP is dictated by: 
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Similarly the time variation of temperature in the 
hadronic phase is given by: 
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In a realistic scenario the value of rj/s may be dif- 
ferent for QGP Hi-lii] and hadronic phases J3JJ37J. 
However, in the present work we take the same value 
of rj/s both for QGP and hadronic matter. 



V. RESULTS 

In case of an ideal fluid, the conservation of en- 
tropy implies that the rapidity density dN/dy is a 
constant of motion for the isentropic expansion [27|j . 
In such circumstances, the experimentally observed 
(final) multiplicity, dN/dy may be related to a com- 
bination of the initial temperature Ti and the initial 
time Ti as Zfrj. Assuming an appropriate value of 
Tj (taken to be ~ 0.6 fm/c in the present case), one 
can estimate Ti. 

For dissipative systems, such an estimate is obvi- 
ously inapplicable. Generation of entropy during the 
evolution invalidates the role of dN/dy as a constant 
of motion. Moreover, the irreversibility arising out of 
dissipative effects implies that estimation of the ini- 
tial temperature from the final rapidity density is 
no longer a trivial task. We can, nevertheless, relate 
the experimental dN/dy to the freeze-out tempera- 
ture, Tf and the freeze-out time, Tf by the relation, 



dN 
dy 



^- = irR 2 A 4a H TJT f /K 
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where Ra is the radius of the colliding nuclei (we con- 
sider AA collision for simplicity) and n is a constant 
~ 3.6 for massless bosons. 

To estimate the initial temperature for the dissi- 
pative fluid we follow the following algorithm. We 
treat Ti as a parameter; for each T^, we let the sys- 
tem evolve forward in time under the condition of 
dissipative fluid dynamics (Eq. [3]) till a given freeze- 
out temperature Tf is reached. Thus Tf is deter- 
mined. We then compute dN/dy at this instant of 
time from eq. U6l and compare it with the experimen- 
tal dN/dy. The value of Tj for which the calculated 
dN/dy matches the experimental number is taken 
to be the value of the initial temperature. Once Tj 
is determined the evolution of the system from the 
initial to the freeze-out stage is determined by the 
Eqs. Q3B EH and QJ 

In Fig. [Tj we display the variation of temperature 
with proper time. It is clear from the results shown 
in the inset (Fig. [TJ that initial temperature for sys- 
tem which evolves with non-zero viscous effects is 
lower compared to the ideal case for a fixed dN/dy. 



Because for a non-viscous isentropic evolution sce- 
nario the multiplicity (measured at the freeze-out 
point) is fixed by the initial entropy. However, for a 
viscous evolution scenario the generation of entropy 
due to dissipative effects contributes to the multi- 
plicity. Therefore, for a given multiplicity (which is 
proportional to the entropy) at the freeze-out point 
one requires lower initial entropy, hence initial tem- 
perature will be lower. It is also seen (Fig. [1]) that 
the cooling of the system is slower for viscous dy- 
namics because of the extra heat generated during 
the evolution. 



A. Photon spectra 

In this section we present the shift in the pt dis- 
tribution of the photons due to viscous effects. The 
integrand in Eq. [3] is a Lorentz scalar, consequently 
the Lorenz transformation of the integrand from the 
laboratory to the co-moving frame of the fluid can 
be effected by just transforming the argument, i.e. 
the energy of the photon (E — pt cosh(y)) in the 
laboratory frame should be replaced by u M p M in the 
co-moving frame of the fluid, where p^ is the four 
momentum of the photon. 

The results presented here are obtained with van- 
ishing bulk viscosity. The effects of viscosity enters 
into the photon spectra through the phase space fac- 
tor as well as through the space time evolution. We 
would like to examine these two effects separately. 
For convenience we define two scenarios: (i) where 
the effects of viscosity on the phase space factor is 
included (6f v is non-zero in Eq. [8]), but the viscous 
effects on the evolution are neglected (rj = in Eq.[5]) 
and scenario (ii) where the effects of r\ ^ are taken 
into account in the phase space factors as well as in 
the evolution dynamics. The space time integrated 
photon yield originating from the QGP in scenario 
(i) is displayed in Fig. [3J Note that the value of 
the initial temperatures for the results displayed in 
Fig. [5] are same (for all rj/s) because the viscous ef- 
fects on the evolution is ignored in scenario (i). The 
viscous effects on the pr distribution of the photons 
is distinctly visible. The higher values of rj/s makes 
the spectra flatter through the px dependence of the 
correction, Sf v . Next we assess the effects of viscos- 
ity on photon spectra for scenario (ii). In Fig. [3] we 
depict the photon spectra for various values of rj/s. 
In this scenario the value of Tj is lower for higher 
rj/s for reasons described above. As a result the en- 
hancement in the photon production due to change 
in phase space factor, 5/^ is partially compensated 
by the reduction in Ti for non-zero rj, which is clearly 
seen in the results displayed in Figs. [5] and [3] 

In Figs.[4]and [5]we exhibit results for the hadronic 
phase for scenarios (i) and (ii) respectively. The ef- 



fects of dissipation on the px distribution of photons 
from hadronic phase is qualitatively similar to the 
QGP phase though the effect is stronger in the QGP 
phase than in the hadronic phase. It is expected that 
the observed shift in the photon spectra due to vis- 
cous effects may be detected in future high precision 
experiments. Finally in Figs. |5] and [7] we plot the 
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FIG. 2: Transverse momentum distribution of thermal 
photons from QGP for various values of rj/s in the sce- 
nario (i). 



Pt spectra of photons for the entire life time of the 
thermal system i.e. the photon yield is obtained 
by summing up contributions from QGP, mixed and 
hadronic phases for different values of rj/s. The ef- 
fect of viscosity for the scenario (i) is stronger than 
than (ii). 
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FIG. 3: Transverse momentum distribution of thermal 
photons from QGP for various values of rj/s in the sce- 
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FIG. 4: Transverse momentum distribution of photons 
from thermal hadrons for various values of rj/s in the 
scenario (i). 



FIG. 6: Transverse momentum distribution of thermal 
photons from the entire evolution history of the system 
for various values of rj/s in the scenario (i). 
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FIG. 5: Transverse momentum distribution of photons 
from thermal hadrons for various values of rj/s in the 
scenario (ii). 



FIG. 7: Transverse momentum distribution of thermal 
photons from the entire evolution history of the system 
for various values of rj/s in the scenario (ii). 



VI. SUMMARY AND DISCUSSIONS 

We have studied the effects of viscosity on the 
evolving QGP produced in nuclear collisions at 
RHIC energies. The generation of entropy due to 
dissipation on the final (experimentally measured) 
multiplicity has been taken into account. The ini- 
tial temperature has been constrained by the mul- 
tiplicity (entropy) at the freeze-out point. The vis- 
cous effects on the photon spectra has been intro- 
duced consistently through the evolution dynamics 
and phase space factors of the particles participating 
in the production processes. It has been noticed that 
the effects of viscosity on the photons from QGP is 
stronger than those originating from hadrons. The 
QGP expected to be formed in Pb+Pb collisions at 



LHC energy will have longer life time and larger vol- 
ume than that of RHIC i. e. the evolution dynamics 
at LHC will be dominated by QGP phase. There- 
fore, the chance of estimating the value of rj/s from 
photon at LHC is brighter. In view of the fact that 
the future experiments are progressing toward preci- 
sion measurement the shift in the pr distribution of 
photons due to dissipative effects at LHC and RHIC 
may be detectable. 

Before closing this section two comments are in or- 
der here. First, as mentioned before for the photon 
production rate from QGP we have used the Comp- 
ton and annihilation processes. We have checked 
that the contribution from these two processes is 
down by a factor of 3-4 compared to the production 
rate obtained from the complete calculation of order 



a s done in Ref. [20j . Taking these higher order pro- 
cesses into consideration in the present scenario in- 
volves a reevaluation of the photon production rates 
with thermal distribution factors containing viscous 
corrections. Secondly, we have confined only to the 
longitudinal flow of the matter in the present work 
ignoring the transverse kick (blue shift) received by 
the photons from radial flow [38|]. However, both 
these factors will affect the photon spectra from ideal 
as well as dissipative scenarios in a similar fashion. 
Therefore, we expect the shift in the transverse mo- 
mentum spectra of thermal photons in the presence 
of dissipative effects which is the main focus of the 
present work, will be detectable even when a more 
rigorous photon production rate along with trans- 
verse expansion is employed 39] . 
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Appendix A: Phase Space 



Considering the energy conservation (£3 = E\ + 
E2 — E) and writing the energies in terms of rapidity 
(Ei = rriiT coshyt) we get, 

E 3 = myr cosh yi + rn-iT cosh y% — px cosh y (22) 

Equating Eqs. |2"T1 and [25] we have, 

m 3T — Vn\ T + m\ T + p T + 2mirm 2 T cosh(j/i — y 2 ) 
-ImyrVT cosh(yi - y) - 2m 2T PT cosh(y 2 - 2/)]f23) 

However, we also have, 



m 3T 



(pl T + ml)2 
= [{Pit + P2T - Pt) 2 + ml}* 

= [PlT + P 2 2T +PT + %PlTP2T COs(0i 2 ) 
-2p T piT COs(4>i) - 2p T P2T cos(0 2 ) 



-m§]* 



(24) 



where, 



In this appendix we derive Eq. 2 from l.The pho- 
ton production rate from the process, 1 + 2 — > 3 + 7 
is given by, 



E 



1 A/" f d 3 pi f d 3 p 2 f d 3 p 3 



dR 



d 3 p 2 (2tt) 8 J 2E 



2E 2 



2E-, 



h(E 1 )f 2 (E 2 )[l±f 3 (E 3 )}\M\^(p 1 +p 2 -p 3 -p()17) 

Performing the d 3 p 3 integration using the delta func- 
tion and using d 3 p/E — prdprdydcj) we get, 



E 



,dR _ 1 AT 
~¥p~ 16 (2tt) 6 



PiTdpiTdyid(j) 1 p2Tdp2Tdy2d^2 



— h(E 1 )f2(E2)[l±f 3 (E 3 )] 



\M\ 2 5(E 1 + E 2 - E 3 - E) 



(18) 



where 0i and 2 are the angles made by the trans- 
verse momenta of first and second particles with 
the transverse momentum of the emitted photon. 
The momentum conservation along the z-direction: 
p 3z = p\ z + p2z — Pz can be written in terms of ra- 
pidity as: 

m 3 T sinh y 3 = myr sinh y\ + rri2T sinh y 2 — Pt sinh £/19) 

Now the energy, E 3 can be written as: 



cos(0i 2 ) = cos(0i) cos(0 2 ) + sin(0i) sin(0 2 ) (25) 

Equating Eq. [23] with Eq.[M] leads to the expression, 

[(Pit cos 0i — pr) cos (f>2 + Pit sin 0i sin 2 ] = 

[(ml + m 2 - m§) + 2mirm2T cosh(j/i — y 2 ) 



2p 



'IT 



-2mi T PT cosh(j/i - y) - 2m 2T p T cosh(y 2 - y) 
-2ptPit cos 4>i] (26) 



Solving Eq. [26] for 2 one gets, 
Pit sin 0i 



tan _i (- 



p\T COS (pi — pr 



) — COS 



H 



2Rp 



21: 



(27) 



where, 

and, 
H 



R = yP\T +Pt~ ZPlTPT COS 01 



(28) 



E 3 = m 3 T cosh 2/3 = \l m\ T + m\ T sinh y 3 (20) 
Substituting Eq.[l9"lin Eq.HOJwe get, 

£3 = ■J [(myr sinh y\ + m,2T sinh y 2 — pr sinh y) 2 + m 3T ] 

(21) " 



(mj + m 2 - m§) + 2mi T m 2 T cosh(yi - y 2 ) 
-2m 1T p T cosh(j/i - y) - 2m 2 TPT cosh(y 2 - y) 
+2prpiT cos 0i (29) 



Now we express the argument of the delta function 
in Eq. fTS] as function of 2 as 

/(0 2 ) = E l +E 2 -E 3 -E 

= miT cosh yi + m,2T cosh y 2 — px cosh y 



- [m 3T + (myr sinh y\ + to 2 t sinh y 2 



-pTsinhy)^ 



(30) 



and performing the 02 integration in Eq. [15] we get, 
,dR 1 M 



E 



d 3 p 16 (2tt) s 



PixdpiT / dp 2 T 



2tt 



dyi / dy 2 

> J-oo JO 

fi(E 1 )f 2 (E 2 )[l±f 3 (E 3 )} 

jMp 

|piTsin(0i - ^2 ) +prsin0 2 Uo 



with the constraint 



H 



'2i?p : 



< 1 



2T 



(32) 



(31) originating from | cos(0) |< 1. 
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